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ABSTRACT 

We use N-body simulations to study the effects that a divergent (i.e. "cuspy") dark 
matter profile introduces on the tidal evolution of dwarf spheroidal galaxies (dSphs) . 
Our models assume cosmologically-motivated initial conditions where dSphs are dark- 
matter dominated systems on eccentric orbits about a host galaxy composed of a dark 
halo and a baryonic disc. We find that the resilience of dSphs to tidal stripping is 
extremely sensitive to the cuspincss of the inner halo profile; whereas dwarfs with a 
cored profile can be easily destroyed by the disc component, those with cusps always 
retain a bound remnant, even after losing more than 99.99% of the original mass. 
For a given halo profile the evolution of the structural parameters as driven by tides 
is controlled solely by the total amount of mass lost. This information is used to 
construct a semi-analytic code that follows the tidal evolution of individual satellites 
as they fall into a more massive host, which allows us to simulate the hierarchical 
build-up of spiral galaxies assuming different halo profiles and disc masses. We find 
that tidal encounters with discs tend to decrease the average mass of satellite galaxies 
at all galactocentric radii. Of all satellites, those accreted before re-ionization (z^6), 
which may be singled out by anomalous metallicity patterns, provide the strongest 
constraints on the inner profile of dark haloes. These galaxies move on orbits that 
penetrate the disc repeatedly and survive to the present day only if haloes have an 
inner density cusp. We show that the size-mass relationship established from Milky 
Way (MW) dwarfs strongly supports the presence of cusps in the majority of these 
systems, as cored models systematically underestimate the masses of the known Ultra- 
Faint dSphs. Our models also indicate that a massive M31 disc may explain why many 
of its dSphs with suitable kinematic data fall below the size-mass relationship derived 
from MW dSphs. We also examine whether our modelling can constrain the mass 
threshold below which star formation is suppressed in dark matter haloes. We find 
that luminous satellites must be accreted with masses above 10 8 -10 9 M Q in order to 
explain the size-mass relation observed in MW dwarfs. 

Key words: galaxies: halos - Galaxy: evolution - Galaxy: formation - Galaxy: 
kinematics and dynamics 



1 INTRODUCTION 

One of the strongest predictions from our present cosmo- 
logical paradigm, Cold Dark Matter (CDM), refers to the 
universal density profile shared by dark matter haloes on 
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all mass scales (Navarro, Frenk & White 1996, 1997; here- 
inafter NFW). These authors suggested a simple formula 
to describe the spherically averaged profile of dark matter 
haloes 



(r/R 3 )[l + (r/R s W 
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This "cuspy" profile diverges towards the centre of the 
halo as p(r) oc r -1 , although the enclosed mass is finite. Sub- 
sequent work has confirmed the absence of central-density 
cores in dark matter simulations (Moore et al. 1999b, Ghigna 
et al. 2000; Fukushige & Makino 2001; Gao et al. 2004; 
Navarro 2004; Diemand et al. 2005; Springel et al. 2008a). 

Despite this consensus, there have been conflicting 
claims about the slope of the inner halo profile, 7 = 
— dlnp/dr. For example, Diemand et al. (2004, 2005) claims 
a central slope of 7 = 1.2, whereas other authors argued 
in favour of steeper, 7 = 1.5, cusps (Moore et al. 1999b; 
Gighna et al. 2000). Recently, Navarro et al. (2010) have 
found small, but significant deviations from an NFW profile 
in the inner-most regions of different Milky Way-size halo re- 
alizations (the Aquarius project; Springel et al. 2008b). This 
result may suggest that the slope of the inner halo profile 
may not strictly be "universal" (see also Gao et al. 2008), 
which might potentially explain the discrepancy found in 
the literature about the value of 7. 

The existence of cores or cusps in the centre of CDM 
haloes has important observational consequences for the dis- 
covery and interpretation of possible signals of dark matter 
annihilation in gamma-ray surveys (Stoehr et al. 2003; Die- 
mand et al. 2007; Kuhlen et al. 2008; Springel et al. 2008b), 
since the strength of the signal goes with the square of the 
dark matter density. 

Interestingly, the presence of cusps in dark matter 
haloes has been challenged by measurements of the rotation 
curves of Low Surface Brightness galaxies (LSBs). Obser- 
vations of LSBs, which are supposed to be embedded in a 
dark matter-dominated potential, seem to favour constant- 
density cores, 7 = 0, for most observed systems (e.g. Flores 
& Primack 1994; Salucci & Burkert 2000; Swaters et al. 2003; 
Simon et al. 2003; Trachternach et al. 2008; Oh et al. 2008; 
de Block 2008, 2010). 

Dwarf Spheroidal galaxies (dSphs) are another exam- 
ple of dark matter-dominated galaxies where significant con- 
straints on 7 can potentially be derived from the kinemat- 
ics of their stellar constituents. Two types of analysis can 
be found in the literature: the first solves the Jeans' equa- 
tions to estimate the mass enclosed within the stellar ra- 
dius. Unfortunately, in pressure supported systems mass and 
orbital anisotropy are degenerate, making any deductions 
of the inner mass profile being cored or cuspy in general 
model-dependent (Koch et al. 2007a,b; Gilmore et al. 2007, 
Evans et al. 2009; Walker et al. 2009). A second analysis 
which aims at breaking the above degeneracy appeals to full 
multi-component distribution-function models, which fit the 
intrinsic and projected moments from hundreds of discrete 
radial velocities. Investigations using this technique seem to 
favour the presence of cores in the dark matter haloes of 
dSphs (Wilkinson et al. 2002; Gilmore et al. 2007). 

The apparent difficulties of CDM to reproduce observa- 
tions on small scales using pure-dark matter simulations has 
not deterred a rich body of work from postulating alternative 
explanations within the present paradigm. Several authors 
have argued that baryons may play a fundamental role in 
sculpting the inner profiles of dark matter haloes, even in 
galaxies that at present are heavily dark matter-dominated. 
For example, mass outflows driven by supernova explosions 
may under plausible assumptions remove dark matter cusps 
(e.g. Navarro, Vincent & Frenk 1996; Read & Gilmore 2005; 



Governato et al. 2010). Also, massive clusters that sink into 
the inner-most regions of dark matter haloes via dynami- 
cal friction will exchange momentum with central particles, 
ejecting them from the cusp and creating a core in the pro- 
cess (Goerdt et al. 2008). 

It is however debatable whether these processes are rel- 
evant on dSph scales, given that only a couple of MW dwarfs 
(Sagittarius and Fornax) hold a bound cluster population, 
and that the large majority of dwarfs with suitable kinemat- 
ics show mass-to-light ratios <; 100. 

Clues to these questions may be gathered from examin- 
ing the overall properties of the satellite galaxy population, 
rather than from individual systems. This paper is aimed 
at studying the influence of dark matter cusps and cores 
on the evolution of satellite galaxies driven by tides, as well 
as exploring whether the satellite population around spirals 
may retain information on the inner structure of dark mat- 
ter haloes. The arrangement of the paper is as follows: Sj2] 
presents the models assumed to describe the host and satel- 
lite galaxies, as well as the numerical set-up. ^presents the 
results of our N-body experiments regarding the influence of 
cusps and cores on the dynamical evolution of dSphs. Ap- 
pendb{X]uses these results to construct a semi-analytic algo- 
rithm that simulates the formation of spiral galaxies through 
the accretion of individual subhaloes. [j4] shows the results 
of applying this code to Milky Way-like galaxies and ana- 
lyzes the properties of the surviving satellite population for 
a range of inner halo slopes and disc masses. We end with a 
brief summary in [J5] 



2 GALAXY MODELS 
2.1 Density profiles 

We assume that dark matter haloes follow a universal den- 
sity profile that can be fitted by the general five-parameter 
formula (Zhao 1996; Kravtsov et al. 1998) 

P{r) = (r/Rs^ll + fr/R^ -<)/- ^ ~ Rvil) ' (2) 

where 7 controls the inner slope of the profile, /3 the outer 
profile, and a the transition between both. In cosmological 
simulations, the characteristic density po and scale radius 
R s are sensitive to the epoch of formation and correlate 
with the halo virial radius via the concentration parameter 

Cvir — Rvir j Rs • 

The virial radius is defined so that the mean over- 
density relative to the critical density is A, 

(4/3)^ lr - Apcrit ' (3) 

where p cr it = 3Hq /StyG, and Ho = 100 h km/s/Mpc is the 
present day value of Hubble's constant. 

The choice of A varies in the literature, with some 
authors using a fixed value, such as NFW, who adopted 
A = 200, and others who choose a value motivated by 
the spherical collapse model, where (for a flat universe) 
A ~ 178fi? n 45 (Lahav et al 1991, Eke et al 1996). The lat- 
ter gives A = 95.4 at z = in the concordance ACDM 
cosmogony, which adopts the following cosmological pa- 
rameters: fi m = 0.3, J1a = 0.7, h — 0.7, consistent with 
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constraints from CMB measurements and galaxy clustering 
(see, e.g., Spergel et al 2007 and references therein). 

The density profile given by eq. [2] may diverge as r — > 
if the inner slope 7 > 0. In the literature such dark matter 
haloes are called "cuspy" , and a prototypical example is the 
NFW profile, which corresponds to (a,/3,-y) — (1,3,1). In 
contrast, "cored" haloes have a converging density profile at 
small radii, i.e. 7 = 0. 

Also, density profiles with outer slopes /3 > 3 lead to 
finite mass models, but for /3 < 3 the cumulative mass pro- 
files diverges as r — > 00. To avoid that problem, Kazantzidis 
et al. (2004) impose a truncation at i? v i r that can be written 
as 



p(r) 



Po 



exp 



Rdc 



where Rd cc is a small quantity, whose value we choose as 

0.1R vlI . 

To obtain a continuous logarithmic slope e is defined as 



-7 - ffgyir 

1 + c« r 



+ 



Rdcc 



(5) 



Our models are purely gravitational and therefore scale- 
free. When convenient for interpretive purposes, we shall 
scale our numerical units to physical parameters assum- 
ing that the host galaxy potential has M v ir = W 12 Mq, 
7?vir = 258 kpc and R s = 21.5 kpc, assuming A = 101. 
These parameters are thus similar to those estimated by 
Klypin et al. (2003) for the MW dark matter halo. 

To explore the effects that the presence of a baryonic 
disc in the host galaxy induces on the satellite galaxy pop- 
ulation, we add a new force component fd = — V$d in the 
form of a Miyamoto-Nagai (1975) disc model 

GM d 



^R 2 + {a + Vz 2 + b 2 ) 2 



(6) 



where Ma is the disc mass, a = 6.5 kpc and b — 0.25 kpc 
are the radial and vertical scale lengths. 



2.2 N-body realizations of satellite galaxies 

Our satellite galaxy models are constructed by generat- 
ing 5-million particle realizations of the profile outlined 
in eq. [2], using the code kindly made available by S. 
Kazantzidis (Kazantzidis 2004, 2006). These models are 
composed purely of dark matter, thus neglecting the con- 
tribution from baryons to the overall potential of satellite 
galaxies. This approximation is valid for the large majority 
of satellite galaxies, i.e. dwarf galaxies, which typically show 
very high mass-to-light ratios (Aaronson 1983, Aaronson & 
Olszewski 1987; Mateo 1998; Walker et al. 2007; Simon & 
Geha 2007; Strigari et al. 2008; Penarrubia et al. 2008a). 

We have chosen parameters for the dwarf consistent 
with those expected for dwarf galaxy halos (see, e.g., 
Penarrubia et al. 2008a). Using the same physical scaling 



adopted for the host galaxy (§ 2.1 1 we consider halo mod- 
els with m v i r = 10 9 M Q , r v - 1T = 25.7 kpc, and c vlI = 23.1. 
In these models the scale radius is independent of the inner 
slope of the profile, r s = 1.11 kpc. We emphasize, however, 
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Figure 1. Density profiles of dark matter halo models for differ- 
ent values of 7, and same virial mass and concentration. Cuspy 
and cored dark matter haloes are those with 7 > and 7 = 0, re- 
spectively. Dotted lines correspond to initial N-body realizations. 
Symbols show the same realizations integrated for 14 Gyr in iso- 
lation. The good agreement between lines and dots implies that 
our numerical setup is adequate and free from numerical artifacts 
induced by numerical limitations down to scales comparable with 
the grid-size of the highest-resolution zone used by SUPERBOX, 
r s /126 ~ 8.8 X 10" 3 kpc. 



that our models are scale free, so that the only relevant pa- 
rameters in the calculations are the relative masses and sizes 
of the host and dwarf. In what follows, capital and small let- 
ters will be used to denote quantities corresponding to the 
"host" or the "satellite" models, respectively. 

Fig.[l]shows the resulting density profiles of the satellite 
for four different values of the slope of the inner density 
profile, 7. Note that differences arise mostly within r < r„. 
Although in those regions the density of dark matter can 
differ by several orders of magnitude depending on 7, in fact 
the differences in enclosed mass are considerably smaller, 
m(< r s )/m vir = 3.86 x 10" 2 ,5.74 x 10" 2 ,8.69 x 10" 2 and 
13.45 x 10" 2 for 7 = 0.0,0.5,1.0 and 1.5, respectively. In 
the following Sections we will see how these relatively small 
differences in central mass play a capital role in the overall 
evolution of satellite galaxies orbiting around spirals. 



2.3 The N-body code 

We follow the evolution of the dSph N-body model in the 
host potential using Superbox, a highly efficient particle- 
mesh gravity code (see Fellhauer et al. 2000 for details). 
Superbox uses a combination of different spatial grids in 
order to enhance the numerical resolution of the calculation 
in the regions of interest. In our case, Superbox uses three 
nested grid zones centered on the highest-density particle 
cell of the dwarf. This center is updated at every time step, 
so that all grids follow the satellite galaxy along its orbit. 
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Each grid has 128 3 cubic cells: (i) the inner grid has a 
spacing of dx = r s /126 ~ 8 x 10 -3 r 3 and is meant to resolve 
the innermost region of the dwarf, where the effects of vary- 
ing 7 are strongest, (ii) The middle grid extends to cover 
the whole dwarf, with spacing r v i r /126. (iii) The outermost 
grid extends out to 50 x r v i r and is meant to follow particles 
that are stripped from the dwarf and that orbit within the 
main galaxy. 

Superbox uses a leap-frog scheme with a constant 
time-step to integrate the equations of motion for each par- 
ticle. We select the time-step according to the criterion of 
Power et al (2003); applied to our dwarf galaxy models, this 
yields At — 4.6 Myr for the physical unit scaling adopted in 
Sec.[2~T1 

2.4 Tests of the dSph N-body model 

We have checked explicitly that our procedure for building 
the dSph model leads to a system that is in equilibrium and 
that does not evolve in isolation away from the prescribed 
configuration. This is shown in Fig. [T] where we plot the ini- 
tial density profiles of four dark matter models with different 
inner slopes (symbols). Dotted lines show the same profiles, 
but after evolving the system in isolation for 14 Gyr with 
Superbox. The good agreement between lines and symbols 
shows that our numerical choices are appropriate, and that 
the evolution of the N-body model, on the scales of interest, 
is free from artifacts induced by the finite spatial and time 
resolution of the calculation. 
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Figure 2. Mass enclosed within the tidal radius (see text) of 
the satellite models presented in Fig. ^ estimated at different 
Galactocentric distances for host galaxies without [upper panel) 
and with a baryonic disc component with a fiducial mass Md = 
0.1M v i r (lower panel), respectively. The "mixed" models refer to 
satellites with cored 7 = 0.0 profiles orbiting about a host galaxy 
with a 7 = 1.0 halo profile. All radii are normalized by the scale 
radius r s of the satellite model. Arrows mark the position of the 
host scale radius, i.e. R s /r s , for this particular satellite model. 



3 TIDAL EVOLUTION OF DARK MATTER 
SUBHALOES 

3.1 Tidal strength 

Satellite galaxies may lose some fraction of their initial mass 
to tides if their orbits bring them close to the host galaxy 
centre. We may quantify the expected strength of the tides 
by computing the "tidal radius" (rtid) of the satellite galaxy 
if this was placed on a circular orbit at pericentre 

(psat(nid)) = 3(pho S t(r pcri )). (7) 

Here (p(r)) denotes mean enclosed density within a ra- 
dius r. Penarrubia et al. (2008b) have shown that the total 
amount of mass, mtid = rn 3 (< rtid), provides an accurate 
estimate of the amount of mass that the satellite may ex- 
pect to retain after completing several orbits in the host 
potential. 

In order to explore the resilience to tides of galaxies 
with different 7- values we have calculated the mass enclosed 
within the tidal radius of the models shown in Fig. [l] as 
a function of pericentric distance using eq. [2] and [7] and 
imposing a universal density profile (but different scaling) 
for both host and satellite models. The results are shown 
in Fig. [2] In absence of a baryonic disc embedded in the 
host halo (i.e. Md = 0), cuspy and cored satellite models 
show similar resilience to tides, although cored satellites are 
slightly more prone to tidal stripping than cuspy ones. Satel- 
lites with divergent profiles show a monotonic decline of m t jd 
as the pericentric radius decreases. Interestingly, this is not 
the case for satellites with a finite central density, which can 



retain a small bound mass fraction m t id ~const even if the 
orbits are fairly eccentric r peT i^r s . 

Motivated by recent claims that star forma- 
tion/feedback processes in the inner regions of dark 
matter haloes can turn cusps into cores (e.g. Governato et 
al. 2010, and references therein) we have also explored a 
"mixed" case wherein our satellite models are cored and 
orbit a host halo with a 7 = 1.0 profile (dotted-dashed 
lines). Fig. [2] shows that in this configuration satellite 
galaxies would undergone tidal disruption if their orbits 
bring them at galactocentric distances 10r s . 

In models where the parent galaxy hosts a baryonic 
disc with a fiducial mass Md = 0.1M v ; r (lower panel) the 
influence of the disc component increases with decreasing 
7 at all pericentric radii. This is because at small orbital 
pericentres, r pcr i J$ R 3 , the mean density of spiral galaxies 
embedded in cored haloes scales as 

(phost) ~ Po + 3aM d /[8Ti(a + b) 3 ]r~^ ri . 

Thus, in models with 7 = 0, the disc component domi- 
nates the host tidal field if r por i 5; (3a[87r(a + fc) 3 ] -1 Md/ po)- 
Complete tidal disruption is assured if the mean density of 
the host is higher than the central density of the satellite. 
This occurs at r pcr j 9aM d /[2(a + b) 3 m s ]r 3 y(c v h-), where 
y(x) = ln(l + x) - x{2 + x)/[2(l + x) 2 ]. Using the model 
parameters adopted in !|2] disruption is expected for orbits 
with r pcr i ^ 32r s . 
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In contrast, at small radii the mean density of spiral 
galaxies with cuspy, 7 = 1, haloes scales as 

(Phost) ~ 3/9 -Rs?-p C 2 ri + 3aM d /[8-K(a + b) 3 }^^. 

Thus in these models the disc tidal field does not dominate at 
small radii. Furthermore, in the limit r por i — > one can read- 
ily show that eq. [7] has an analytical solution, which can 
be written as r t = [m s i? s /(c v i r )/(3r s M v ir/(Cvir)] 1/ ' 2 r pcr i, 
where f(x) = ln(l + x) — x/(l + x). Since at small radii 
m s (< r t id) oc r t id, this implies that m s (< r t i<j) = if and 
only if r per i = 0, which suggests that only radial orbits may 
lead to full tidal disruption of satellites orbiting in such 
a galaxy potential. A similar prediction was made by Go- 
erdt et al. (2007) using different analytical arguments. As 
we shall see below, the apparent "indestructibility" of cuspy 
subhaloes is confirmed by our N-body simulations. 

It is clear from inspecting this Figure that baryonic discs 
may represent an important factor in the survival of satellite 
galaxies around spiral galaxies, and that the disc influence 
strongly varies depending on the slope of the inner density 
profile. 



3.2 The orbits 

To explore self-consistently the tidal evolution of satellite 
galaxies we place the N-body models outlined in §2.2| in or- 
bit within the host galaxy assuming different combinations 
of apocentric/pericentric ratio. All orbits have apocentres 
r apo ~ 0.70i? v ir = 180 kpc. We consider a large range of 
pericentric radii r- pcr i = 1/2,1/5,1/20 and 1/100 of r apo . 
Note that the smallest pericentre is comparable to the scale 
radius of the satellite, and that this type of close-range en- 
counter is rare in self-consistent cosmological simulations. 
In this contribution, this series of orbits can be regarded as 
one of increasing tidal strength for eccentric orbits of simi- 
lar orbital period. In all cases we follow the evolution of the 
satellite for tn = 14 Gyr in the scale units introduced in 
§2.1[ which corresponds to 3-5 orbital periods. 



3.3 Mass loss 

Tides are strongest at pericentre and trigger mass loss 
episodes that recur with every pericentric passage. The pro- 
gressive stripping of satellites is shown in Fig. [3j where dots 
represent the fraction of mass that remains bound to the 
satellite as a function of time. Pericentric passages are eas- 
ily recognizable as the times when sudden drops in bound 
mass occur. Between pericentric passages the bound mass 
fraction remains nearly constant, suggesting, as indicated 
above, that tidal forces act in an impulsive regime and have 
little influence on the main body of the satellite at times 
other than pericentre. 

As expected from Fig. |2j in the absence of a disc com- 
ponent in the host galaxy (left panel) the difference in mass 
loss that results from varying the slope of the inner density 
profile of our satellite models is in general small, and mostly 
arises at the latest stages of the orbital evolution. An excep- 
tion can be gleaned from the most-eccentric modes plotted 
in Fig. [3] (r pcr i : r apo = 1 : 100), which show that after 
the first pericentric interaction the bound mass fraction of 
satellites with prominent cusps [filled circles) is a factor 2-3 



lower than that of the cored models [open squares). However, 
this situation reverses as the number of pericentric passages 
increases, which means that, when cored haloes begin to 
disrupt, they do so more rapidly than cuspy haloes. 

The right panel of Fig. [3] shows that the presence of a 
disc component in the host introduces remarkable differences 
between the mass evolution of cored and cuspy satellites. 
Whereas the rate at which satellites with prominent cusps, 
7 = 1.5, lose mass is barely altered, the full tidal disruption 
of models with a 7 = 0.0 profile is assured if their orbits 
bring them into the vicinity of the disc. 

Indeed, dark matter haloes with 7 > 1.0 show an ex- 
treme resilience to tides in all models explored here. As Fig. [3] 
shows, the first pericentric passage results in a significant de- 
crease of the bound mass fraction. However, the mass loss 
rate progressively slows down as the number of consecutive 
pericentric passages increases, and the bound mass tends 
to converge to a finite value. This result implies that dwarf 
galaxies embedded in haloes with strong cusps, 7 > 1.0, can- 
not be fully disrupted by the tidal field of the host galax y, as 
expected from the analytical arguments given in pin This 



has interesting consequences for the properties of the satel- 
lite galaxy population, an issue to which we return below. 

The tidal mass stripping of dark matter haloes pro- 
ceeds gradually from the outside in (e.g. Stoehr et al. 
2002; Kazantzidis et al. 2004; Kampakoglou & Benson 2007; 
Penarrubia et al. 2008b) . This is illustrated in Fig. |4j where 
we plot, at consecutive apocentres of the orbit, the position 
of the satellite particles in the satellite frame. This clearly 
shows how a halo is gradually stripped from the outside in, 
following an "onion-peel" process in which mass loss affects 
predominantly the outer, low-density regions of the satellite 
galaxy. 

By colour-coding the particles according to the value 
of the projected surface density at their location we can 
appreciate fundamental differences in the evolution of dark 
matter haloes depending on whether they are cuspy (left 
panel) or cored (right panel), even if both satellite models 
have lost the same fraction of their original mass, 97%, after 
three pericentric passages. The first interesting difference is 
that the spatial sizes of the cuspy satellite models suffer a 
dramatic decrease as a result of mass loss. Yet the bound 
remnants show densities that have been barely altered by 
the action of tides. The response of cored satellite models to 
tidal mass loss is practically the opposite. In these systems 
the change in size if much less pronounced, and the density 
drops more strongly than in the cuspy models at all radii, an 
issue that we explore in more detail in the following Section. 



3.4 Evolution of the density profile 

Fig.[5]shows the evolution of the density profile of cored (7 = 
0.0, open circles) and cuspy (7 = 1.0, filled circles) satellite 
models for the most extreme orbit (r por i : r apo = 1 : 100) 



1 Note that Hayashi et al. (2003) stated that, if haloes are trun- 
cated at a radius r < 0.77r s , that would result in the total dis- 
ruption of the remnant due to the positive binding energy of the 
NFW in that region. However, in |3.4| we show that tidal stripping 
does not truncate the density profile of dark matter haloes (see 
also Penarrubia et al. 2008b, 2009) 
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Figure 3. Bound mass fraction as a function of time for different orbits and values of 7. Left panels show the mass evolution of subhaloes 
in a host halo without disc. Right panels show the effects of injecting a disc component in the host galaxy potential. Note that the disc 
influence on the survival of dark matter subhaloes strongly depends to the slope of the inner density profile, 7. In particular, the mass 
evolution of satellite galaxies with prominent cusps, i.e. 7 = 1.5, is barely altered by the disc presence, whereas cored satellites, i.e. 
7 = 0.0, suffer severe mass loss if their orbits bring them into the disc vicinity. 



probed in our series. The systems are shown at three consec- 
utive orbital apocentres, when the satellites have had some 
time to find dynamical equilibrium after being perturbed by 
the host tidal field at pericentre. The dotted/dashed lines are 
the same in each panel and show the initial density profiles 
for ease of reference. 

A few interesting points can be gleaned from this Fig- 
ure. One is that it is possible to identify two distinct regions 
on the basis of the shape of the density profile: an inner 
region where the profile varies smoothly with radius, and 
an outer region where an "excess" of material is present 
with respect to a naive extrapolation of the inner profile. 
Several authors (e.g Oh, Lin & Aarseth 1995; Mayer et al. 
2002; Johnston et al. 2002) have shown that this type of 
"breaks" in the density profile are transient features that 
arise from the presence of unbound particles in the out- 



skirts of the system as they escape from the satellite. More 
recently, Penarrubia et al. (2009) have shown that the po- 
sition of the "break", rb, marks the radius at which the lo- 
cal crossing time exceeds the time elapsed since pericentre, 
t — t p , so that rb oc a(t — t p ), where a is the mean velocity 
dispersion of the satellite. Particles inside the break radius 
are thus within the region where enough time has elapsed 
for equilibrium to be reached. 

A second interesting point of note is that the evolution 
of the density profile is remarkably sensitive to the shape of 
the inner profile. Tides acting on cuspy haloes tend to induce 
strong drops of density in the outer regions, as expected from 
Fig. [4] whilst keeping the shape of the inner profile barely 
altered, even after the satellite has been stripped of 97% of 
its initial mass. The bottom-right panel shows that, for a 
similar mass loss, tides tend to lower the central density of 
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Figure 4. Projection onto the orbital plane of a cuspy (7 = 1.0, left panels) and cored (7 = 0.0, right panels) satellite models that 
follow a highly eccentric orbit (r por ; : r apo = 1 : 100) about a host galaxy with no disc component. Each of the panels correspond to 
consecutive orbital apocentres for each model. After three consecutive pericentric passages, both satellites have lost 97% of their initial 
mass. Particles have been colour-coded according to the value of the projected surface density at their location. Note how in both models 
the stripping proceeds from the outside in, and that the cuspy model shrinks at a faster pace than the cored one, without a substantial 
drop of the central density. 



cored haloes more severely than in the case of haloes with a 
centrally-concentrated profile. 

The evolution of the satellite's circular velocity curve 
has a different amplitude depending on the slope of the inner 
halo profile. This is shown in Fig. [5j in which we mark the 
position where the circular velocity reaches its peak,r max , 
at consecutive orbital apocentres . As this Figures shows, 
mass loss results in a progressive shift of r max towards the 
satellite centre, a process that is more pronounced in the 
case of cuspy satellite models. We inspect this issue in some 
detail in the following Section. 

3.5 Evolution of structural parameters 

As shown by Hayashi et al. (2003) and Pefiarrubia et al. 
(2008b), the evolution of the structural parameters of NFW 
haloes as they are tidally stripped does solely depend on the 
total amount of mass lost, and not on the details of how the 
mass is stripped. 

Fig. [6] shows that this result is generic to dark mat- 
ter haloes that follow the density profile outlined by eq. 
The left and right panels of the Figure show respectively 
the fractional change of the peak velocity, v m&K , and its lo- 
cation, r max , as a function the satellite bound mass fraction 
/bound = m s /m s (t = 0) for three of the four orbits surveyed 
in our study. Different symbols correspond to different orbits 
and are shown at consecutive orbital apocentres when the 
main body of the satellite galaxy is close to dynamical equi- 
librium. Open and filled symbols denote models in which a 
disc component has been embedded in the parent galaxy. 

This Figure illustrates a few interesting differences in 
the evolution of satellite galaxies depending on how dis- 
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Figure 5. Halo density profile of the cuspy and cored halo mod- 
els shown in Fig. [4] For comparison we also plot the initial den- 
sity profiles in each panel, where dotted and dashed lines denote 
haloes with 7 = 1.0 and 0.0, respectively. The quantity /bound 
refers to the bound mass fraction. Arrows indicate the position 
of rmax, i-e. the radius where the circular velocity curve peaks. 
Note that the density profile in the inner regions of dark matter 
haloes is preserved even when large fractions of the initial mass 
are lost to tides. 
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tended the halo profile is in the inner-most regions. The 
first is that, at a given mass loss fraction, the peak veloc- 
ity declines more significantly as we decrease the value of 
7. The second difference is in the evolution of the satellite 
size, which is clearly more substantial for satellite galaxies 
with cusps, as expected from Fig. [4] and [5] The right panel 
of Fig. [5] shows that r max is severely shortened by the ac- 
tion of tides if the slope of the inner halo profile is 7 > 0. 
Interestingly, the location of the peak velocity of cored mod- 
els cannot drop further than a factor ~ 3, regardless of the 
amount of mass lost to tides. This is in marked contrast 
with cuspy haloes, which may shrink into oblivion if tidal 
stripping carries on indefinitely. 

Irrespective of pericentre, number of completed orbits, 
host potential, or slope of the inner density profile, the evo- 
lution of the structural parameters depends solely on how 
much the satellite mass has varied. This is a remarkable re- 
sult, especially because our series of models includes extreme 
cases in which the satellite models lose more than 99.99% 
of their original mass. Following Penarrubia et al. (2008b), 
we use a simple empirical formula to fit the evolution of the 
structural parameters 



9{x) = 



2V 



(8) 



where x = m s /m 3 (t — 0) and g(x) represents either u max 
or r max . The best-fit values of \i and r\ have been plotted in 
each of the panels of Fig. [6] and the results are shown with 
dashed lines. 

The fact that the evolution of satellite galaxies as they 
are stripped by tides follows well-defined tracks will be used 
in the following Sections to build a simple semi-analytical 
algorithm that follows the disruption of dark matter haloes 
with different inner profiles with reasonable accuracy (see 
Appendix A). 



4 SEMI-ANALYTIC REALIZATIONS OF 
GALACTIC MERGER TREES 

The results presented in ^3] show that the response of galax- 
ies to tides strongly depends on the slope of the inner density 
profile of dark matter haloes, 7, as well as to the presence 
of a baryonic disc embedded in the host galaxy. This raises 
the question of whether the satellite populations surround- 
ing spiral galaxies may provide insights into the inner struc- 
ture of dark matter haloes. In order to explore this issue, 
simulations that follow the hierarchical formation of spiral 
galaxies through the accretion of individual satellites are re- 
quired. We begin by outlining a method to construct merger 
trees of spiral galaxies, and then discuss the effects of vary- 
ing 7 and the disc-to-parent halo mass ratio on the present 
population of satellite galaxies. 



4.1 Build-up of merger trees 

4-1.1 The host galaxy 

We construct merger trees using Monte Carlo methods. 
Specifically we employ the merger tree algorithm described 
by Parkinson et al. (2008) which is itself an empirical modi- 
fication of that described by Cole et al. (2000). We adopt the 
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Figure 6. Evolution of the peak circular velocity f m ax = 
V c (r max ) and its location, measured at the apocentres of the var- 
ious orbits in our simulation series. Open and closed symbols 
denote simulations with and without a baryonic disc component 
(Mj = 0.1M v i r ), respectively. The dashed lines are our fits to the 
evolutionary tracks using eq. [8]. 
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Figure 7. The evolution of the host galaxy parameters for five 
different merger tree realizations of a halo with Af v ; r = 10 12 AfQ 
at z = is shown in the upper panels. Here we plot the virial mass 
and concentration in the upper-left and -right panels, respectively. 
The disc component embedded in the host galaxy also grows with 
time. For this particular plot we show the evolution of a disc with 
Af,j = 0.1Af v i r and a = 6.5 kpc at z = 0. The disc mass and size 
are plotted in the lower-left and -right panels, respectively. 
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parameters (Go, 71, 72) = (0.57,0.38,-0.01) that are found 
to provide the best fhrl to the statistics of halo progenitor 
masses measured from the Millennium Simulation by Cole 
et al. (2008). We use a mass resolution (i.e. the lowest-mass 
halo which we trace in our trees) of W^h -1 Mq, which is 
sufficient to achieve resolved galaxy properties for all of the 
calculations considered in this work. 

Within this evolving hierarchy of dark matter, the evo- 
lution of the host galaxy is followed using the Galform 
semi-analytic model of galaxy formation (Cole et al. 2000). 
In particular, we use the implementation of Bower et al. 
(2006) which includes feedback from active galactic nuclei 
and which produces a good match to the evolving stellar 
mass function of galaxies. For the purposes of this work the 
key results extracted from the model are the evolving prop- 
erties (mass and size) of the disc forming at the centre of 
the dark matter halo. These are computed in Galform by 
determining, at each time step, the rate at which hot diffuse 
gas in the halo is able to cool and sink towards the halo cen- 
tre, where it becomes rotationally supported, forming a disc. 
Simple prescriptions for star formation and feedback from 
supernovae then allow the gaseous and stellar mass of the 
galaxy to be determined. Knowledge of the angular momen- 
tum of the disc, assumed to have been conserved throughout 
the process of cooling, collapse and disc formation, allows for 
an estimate of the size to be found, following the methods 
of Mo, Mao & White (1998). 

We run ten different realizations of a MW-like galaxy 
with a mass of 10 12 Mq at z = 0. Although all models end up 
with the same mass and virial concentration, the assembly 
of these haloes is remarkably distinct, a direct manifestation 
of cosmic variance. A useful distinction between different 
halo models can be made by deriving the time at which a 
significant mass fraction has been accreted. For example, 
Fig. [7] includes examples of late-forming galaxies, which do 
not accrete 10% of their present mass until redshit z ~ 2 
(e.g solid lines), as well as early-forming galaxies, which have 
this mass fraction in place by z ~ 5 (e.g. long-dashed and 
dotted-dashed lines). 

The disc mass is set by choosing a value for the present 
disc-to-halo mass ratio £ = (^Md/ . Since the bary- 

onic content of galaxies cannot exceed the universal baryon 
fraction, fi;,/fi m ~ 0.17 (Spergel et al. 2007), we consider 
disc masses within the range < £ < 0.2. At any given red- 
shift, the disc mass is re-scaled as M' d (z) — Md(z)/Md{z = 
0)£Af v i r . In general the disc component forms at z^2 (e.g. 
solid lines), although in some extreme cases the formation 
process can start as early as z ~ 5, (e.g. long-dashed lines). 
The disc mass suffers strong fluctuations with time, as shown 
in the bottom-left panel of Fig. [7] The formation time of 
discs plays a key role in the mass evolution of satellites, as 
discussed in the following Sections. 



2 Benson (2008) found an alternative set of parameters which 
provided a better match to the evolution of the overall halo mass 
function but performed slightly less well (although still quite well) 
for the progenitor halo mass functions. We have chosen to use 
the parameters of Parkinson et al. (2008) as we wish to get the 
progenitor masses as correct as possible. 
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Figure 8. Cumulative number of satellite galaxies accreted with 
fis(2 aC c) > 1O S M0 as a function of redshift {upper panel) ex- 
tracted from the five merger-tree realizations shown in Fig. [7] 
Distance from the host galaxy centre at the time of accretion, 
r acc , as a function of redshift (lower panel). Note that if orbital 
decay did not take place, r acc ~ r apo throughout the orbital evo- 
lution of satellite galaxies. 

4-1-2 The satellite population 

Our merger tree realizations follow the evolution only of 
subhaloes with masses above 1O 8 M0 when they first cross 
the host virial radius. This choice is motivated by the fact 
that primordial gas is unable to cool by Lyman-Q emis- 
sion and initiate star formation in haloes with a temper- 
ature below 10 4 K, which approximately corresponds to the 
above mass threshold. It is typically assumed that lighter 
subhaloes remain dark (e.g. Haiman et al. 2000), and thus 
undetectable through photometric surveys. Unfortunately, 
determining the number of stars that a given subhalo can 
form from primordial gas is still a topic of hot discussion 
(see Tassis et al. 2008; Ricotti 2009 and references therein). 

Our calculations do not include a baryonic component 
in the accreted subhaloes and may therefore be regarded as 
pure dark matter models. Despite the uncertainties in the 
properties of the stellar components embedded in our satel- 
lite models, the following Sections will show that a number 
of general inferences may be made in reference to the struc- 
ture of satellite galaxies. 

The cumulative number of satellite galaxies accreted at 
a given redshift, 7V acc , is shown in the upper panel Fig. [8] 
Since the host halo grows by accreting dark matter sub- 
haloes, the number of accreted satellites and the virial mass 
are closely correlated at any given epoch. An interesting 
point of note is that a host galaxy with M v i T (z = 0) = 10 12 
is expected to accrete around ~ 10 3 satellites that may po- 
tentially host a stellar component. This is in remarkable con- 
trast with the number of known satellite galaxies in the MW, 
which to date adds up to a mere couple of dozens (e.g. Be- 
lokurov et al. 2007; Koposov et al. 2009). The apparent mis- 
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match between theoretical predictions and observations is 
typically referred as the "missing satellite problem" (Moore 
et al. 1999a; Klypin et al. 19990 

When a subhalo first crosses the host virial radius, it be- 
comes a satellite galaxy. The inner satellite structure is fixed 
by applying the virial mass-concentration relation found by 
NFW. Its orbital parameters are drawn from the distribu- 
tion of Benson (2005), which was measured from N-body 
simulations. This distribution gives the radial and tangen- 
tial velocity components of the orbit. By definition, the dis- 
tance from the host galaxy centre at accretion time, r acc , 
equals the host virial radius, and closely correlates with or- 
bital apocentre (Ludlow et al. 2009). Thus, a natural result 
of the inside-out growth of galaxies is a strong correlation 
between the accretion time of satellites and their location 
within the host. 

Indeed, satellites accreted at early times can only be 
found in the inner-most regions of their parent galaxies. This 
is illustrated in the lower panel Fig. |8j where we show the 
infall distance as a function of redshift for the merger-tree 
realizations shown in Fig. [7] Two distinct epochs can be dis- 
tinguished on the basis of the effects introduced by cosmic 
variance: a late epoch, i.e. during the last 8 Gyr (z 1), when 
the orbits of accreting satellites have a similar distribution 
of orbital parameters, and an early epoch (z^, 1), when the 
accretion history of the host galaxy is dominated by cosmic 
variance. This result suggests that the properties of the sur- 
viving satellite galaxy population at z — will show some 
scatter between different merger trees, and that that scat- 
ter will be more evident for satellites that were accreted at 
early, z£l, epochs. 

4.2 Spatial distribution of subhaloes 

The destructive action of tides is illustrated in Fig. [9] where 
we show the present mass and galactocentric distance of 
cuspy (left panels) and cored (right panels) satellite models 
extracted from the ten host galaxy realizations outlined in 
§4.1| For ease of reference, we mark with dashed lines the 
minimum mass that satellites must carry as they are ac- 
creted into the host in order to be included in our merger 
trees. Thus in this Figure all satellites with m s < 10 8 Mq 
have lost some fraction of their original mass to tides. 

Fig. [9] illustrates a number of interesting points. The 
most obvious is that satellites accreted at early times tend 
to lose a large fraction of their original mass to tides, which is 
highlighted in this plot by marking the masses and positions 
of galaxies that fell into the host galaxy at z > 6 (solid dots). 
This redshift roughly corresponds to the epoch at which the 
Universe was completely re-ionized (e.g. Fan et al. 2002), 
and thus we refer to these galaxies as re-ionization fossilx~\ 
It is typically believed that such systems went through a 
short phase of star formation before they crossed the virial 

3 Recent contributions show that the discrepancy between the 
predicted and observed number of satellites in the Milky Way 
alleviates after taking into account the detection biases as well 
as volume corrections of SDSS data (e.g. Tollerud et al. 2008; 
Koposov et al. 2008). 

4 Note that in Ricotti (2009) this term has a different meaning: 
in particular, it refers to dwarf galaxies with peak velocities below 
20 km/s throughout their whole orbital evolution 



radius of the parent galaxy and lost their primordial gaseous 
component through ram pressure stripping, a process cap- 
tured by the N-body simulations of Mayer et al. (2007). The 
stellar compositions of these systems are thus expected to 
be distinct from those of the satellite galaxies that contin- 
ued forming stars after re-ionization and were accreted at 
a later time. Among the known MW satellites Belokurov et 
al. (2009) identifies four possible candidates for being fossils 
of re-ionization: Segue I and II, Bootes II and Coma. Future 
spectroscopic surveys of MW dSphs will possibly extend this 
census. 

Because the initial orbital apocentres of our satellite 
models roughly correspond to the virial radius at accretion 
time, which scales as R V1T tx (1 + z) , re-ionization fos- 
sils are expected to populate the inner-most regions of the 
host, r<R vil (z = 0)/7<30 kpc (see also Gao et al. 2009). 
Whether or not re-ionization fossils survive tidal disruption 
to the present day strongly depends on the inner profile of 
dark matter haloes, as well as on the mass of the host disc. 

Fig. [9] shows that the disc component strongly depletes 
the inner regions of the host galaxy from satellites. Due 
to their small apocentres, re-ionization fossils count among 
the satellites most prone to undergo full tidal disruption in 
7 = 0.0 models. As a result, the number of surviving re- 
ionization fossils is close to negligible if the inner profile of 
dark matter haloes is cored. Although cuspy satellite mod- 
els also shed large fractions of their mass to tides, a bound 
remnant always survives the shredding process (see Fig. [3}. 
In these models the main effect of the disc component is to 
lower the present mass (and luminosity) of these systems. 

A second point of interest is that the satellite galaxies 
that have have been acted on by tides extend all the way 
out to the host virial radius. This is because in CDM most 
satellites move on highly eccentric orbits (e.g. Benson 2005). 
Thus, a significant number of galaxies with relatively large 
apocentres are likely to reach the inner-most regions of the 
host and experience severe stripping during a Hubble time, 
as shown in Fig. [9] 

In Fig. [To] we show the mean cumulative number of 
satellites per realization (upper panels) and bound mass 
(lower panels) as a function of galactocentric distance, as 
derived from the above models. Thin and thick lines distin- 
guish between fossils of re-ionization and the whole satel- 
lite sample, respectively. As anticipated above, the cumula- 
tive number of cuspy satellites barely depends on disc mass, 
again illustrating the fact that cusps cannot be completely 
destroyed by tides. Cored satellites, in contrast, can be fully 
disrupted by tides, which translates into a progressive de- 
crease of the total number of satellites as more massive discs 
are assumed. This is particularly obvious for the satellites ac- 
creted prior to re-ionization (i.e. z acc <;6). Adopting a fidu- 
cial disc mass Md = 0.1M v ; r and a cored halo profile, we 
expect of the order 2-3 fossils of re-ionization in a MW-like 
galaxy. This number nearly quadruples in halo models with 
7 = 1.0 and is barely dependent on disc mass. 

A result that may have important observational conse- 
quences for the detection of dark matter annihilation signals 
in gamma-ray surveys refers to the large abundance of satel- 
lite galaxies at small galactocentric radii. Cuspy models pre- 
dict the presence of the order of 100 satellites within 10 kpc 
from the host centre with masses that range between 10 3 - 
1O 8 M0, a number that barely depends on the disc mass. 
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Figure 9. Satellite mass versus Galactocentric radius at z = for satellite models with m s (z acc ) > 1O 8 M0 and different dark matter 
profiles. Here we stack all the models from ten different realizations of a host galaxy with M v - lT = 1O 12 M0 and r v j r = 210 kpc at z = 0. 
Note that all subhaloes had masses above the dashed lines at their time of accretion. Closed symbols denote satellites that were accreted 
into the host halo with virial masses above 1O 8 M0 at z > 6 (the so-called fossils of re-ionization) . Note that injecting a baryonic disc in 
our host galaxy models has the effect of depleting the number of fossils of re-ionization if haloes are cored (7 = 0.0), whereas if haloes 
are cuspy (7 = 1.0) the main effect is an increase in the mass scatter at a fixed galactocentric radius. 



This number drops by a factor ~ 10 in models that as- 
sume a cored dark matter profile and a fiducial disc mass 
of Md = O.lMvir- Since the gamma-ray flux from a particu- 
lar source is inversely proportional to the distance squared, 
these galaxies may provide interesting targets for dark mat- 
ter annihilation surveys, an issue that we will explore in a 
separate contribution. 

The lower panels of Fig. [lO] show that the variation of 
the mean satellite bound mass as a function of galactocen- 
tric distance could potentially constrain the inner profile of 
satellite galaxies. In particular, our models show that the 
satellite mass is expected to drop at distances r R s ~ 25 
kpc if haloes have a cuspy profile. This is in a marked con- 
trast with cored models, where the effect is the opposite. 
Here the negative mass gradient arises simply from the de- 



pletion of satellites with small orbital apocentres. Satellites 
that at z = are located at distances comparable to the 
disc size happen to be observed by chance near pericentre. 

An interesting point concerns the extreme stripping 
suffered by a large number of satellites in these models. 
Penarrubia et al. (2008a,b) estimate that satellite galaxies 
must lose ^ 90% of the dark matter envelope before stars 
begin to be shed to tides. In our models this implies that all 
satellites in Fig. [9] with masses ^ 10 7 Mq have lost a frac- 
tion of their stellar components to tides and have associated 
stellar streams. The fact that tidal stripping is very sensi- 
tive to the inner profile of dark matter haloes suggests that 
surveys of stellar streams in spiral galaxies may put strong 
constraints on the inner profile of dark matter haloes. Unfor- 
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Figure 10. Mean cumulative number of satellite galaxies (up- 
per panels) and mean bound mass (lower panels) as a function 
of distance from the host centre obtained from the models shown 
in Fig. [9] We measure these quantities for two satellite samples: 
those that were accreted at z > 6, to which we refer as fos- 
sils of re-ionization (thin lines), and the whole sample, which in- 
cludes all satellites with masses above 10 s Mq at accretion (thick 
lines). Solid, dotted and dashed lines denote host galaxies with 
Md/M v j r = 0.0,0.1 and 0.2, respectively. Note the present mass 
and spatial distribution of satellite galaxies strongly depends on 
the slope of the inner halo profile. 

tunately, the absence of a stellar component in our satellite 
models prevents us from exploring this issue in depth. 

4.3 The masses of Milky Way and M31 dSphs 

The results shown in previous Sections suggest that by sim- 
ply counting the number of satellites around spiral galaxies 
one could make strong inferences on the inner structure of 
CDM haloes. However, this requires theoretical models to be 
able to predict properties of the stellar populations hosted 
by dark matter subhaloes and to make comparisons that in- 
volve observable quantities, such as galaxy luminosities and 
surface brightnesses. In practice, this is difficult because, in 
spite of strong theoretical efforts (e.g. Maccio et al. 2009; 
Cooper et al. 2009), the processes that shape the baryonic 
properties of satellite galaxies during their formation and 
later accretion are still rather uncertain. 

A more promising way forward might arguably be to 
compare the masses of satellites derived from the kinemat- 
ics of the stellar tracers against theoretical predictions. In 
MW-size galaxies, the vast majority of satellite galaxies cor- 
respond to dwarf spheroidals (dSphs). Although dSph mass 
estimates are often dominated by large uncertainties (e.g. 
Kravtsov 2010 and references therein), recent contributions 
have shown that the satellite mass enclosed within both the 
projected and tri-dimensional half-light radii of these sys- 
tems is well constrained and robust to a wide range of halo 
models and velocity anisotropies (Walker et al. 2009, Wolf et 



al. 2009) . These data exist for most of the known MW satel- 
lites and a handful of M31 ones, which altogether encompass 
half-light radii that span nearly two orders of magnitude in 
size. 

In Fig. [II] we show the mass estimates for MW (Walker 
et al. 2009) and M31 (Collins et al. 2009; Kalirai et al. 2009) 
dwarfs with closed and open symbols, respectively. Lines 
are derived from the sample of models shown in Fig. [9] by 
simply measuring the mass of individual satellites enclosed 
within the range of radii spanned by the observational data, 
ms(< r). Error bars attached to the curves denote devia- 
tions around the mean, whereas those associated with the 
symbols represent observational uncertainties. Some of the 
surviving satellite galaxies in our models have smaller sizes 
than the largest half-light radii in the sample. To remove 
this bias, only haloes with r v - lr > 2r are considered in the 
calculation of m s (< r). 

This Figure shows two remarkable results. The first one 
is that cuspy satellite models (7 = 1.0) provide a better de- 
scription of the masses of Local Group dwarfs than cuspy 
(7 = 0.0) ones. This is particularly obvious at small radii, 
i.e. the region populated by ultra-faint dSphs, where the dif- 
ference between both halo models are largest (see Fig. [I]). 
Although not shown here, we have also explicitely confirmed 
that this mismatch do not arise from biases between the 
satellite masses and their galactocentric distance or their 
time of accretion. It is also worth noticing that the devia- 
tions around the mean are relatively small taking into ac- 
cout the large range of bound masses encompassed by our 
satellite models, and provide a good representation of the 
scatter shown by the observational data in models where 
the disc component represents a small mass fraction of the 
host halo. This result naturally accounts for the appearance 
of a "universal mass profile" for dwarf galaxies (Walker et 
al. 2009). 

The second interesting point concerns the impact of the 
disc component on the masses of satellite galaxies. Clearly, 
the mean enclosed mass at a given radius tends to decrease 
as more massive disc models are adopted. At a given ra- 
dius, the mean satellite mass can drop as much as a factor 
~ 3 for the most massive disc components considered here, 
Md = 0.2M v i r . Interestingly, the systematically low masses 
shown by the handful of M31 satellites with available kine- 
matics may be the telltale evidence of a relatively heavy 
disc component. In this scenario, the MW and M31 may be 
embedded in similar-size haloes, but the disc-to-halo mass 
ratio of M31 would be a factor ~ 2 higher than that of the 
MW. 

Interestingly, by converting K-band absolute magni- 
tudes into stellar masses Hammer et al. (2007) have recently 
estimated that the disc of M31 may approximately be a fac- 
tor ss 2 more massive than that of the MW, although these 
measurement still have large uncertainties. This result, if 
confirmed, would lend support to the scenario that explains 
the relatively low mass of M31 satellites as tidal in origin. 

Further clues to the mass of M31 disc may be gathered 
from the faint-end of the dSph population. Our modeling in- 
dicates that the mass difference between M31 and MW dSph 
populations should also be present in dwarfs with half-light 
radii smaller than those shown here. This CDM-motivated 
prediction may be verified with present observational capa- 
bilities. 
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Figure 11. Mass enclosed within fixed radii in the satellite frame 
for the models shown in Fig. [9] Lines denote average values and 
error bars deviations around the mean. Filled and open symbols 
show, respectively, the masses of the MW and M31 dwarf galaxies 
with known velocity dispersion estimated at the half-light radius, 
^halfi using Walker et al. (2009) formula. To allow a meaningful 
comparison against observational constraints, only satellite mod- 
els with r v ; r > 2r are considered in the calculation (see text). 



4.4 Mass Threshold for star formation 

In previous Sections a mass threshold at m t hrea = 
wi s (2acc) min > 10 8 was imposed to account for the fact 
that only satellites with masses m s > m t hres can form stars 
(e.g. Heinman et al. 2000). There is, however, some contro- 
versy about the value of m t hres that defines the minimum 
galaxy mass (see e.g. Tassis et al. 2008 and Ricotti 2009 
for disparate results that place m t hr C s above and below the 
adopted 1O 8 M0, respectively). 

Because the satellite mass function diverges as m s — > 0, 
imposing a mass cut onto the satellite population implies 
that satellites with m 3 ~ m t hros dominate in number. As a 
result, increasing (decreasing) the value of the mass thresh- 



11 



old shifts upward (downward) the curves shown in Fig 
at all radii. This suggests that the masses of Local Group 
dwarf galaxies can effectively be used to place constraints 



on m thr . 



Fig. |12| shows a crude attempt to estimate the minimum 
satellite mass above which star formation is expected to oc- 
cur. To estimate the accuracy of our models in describing 
the masses of the MW dwarf galaxies shown in Fig. |11| we 
apply ax 2 fit, where 



E 



N obs 



E 



1 (Vo 



(9) 



and y = log 10 M. The observational variance is a b s = 
(?m/{M log 10 ), where M denotes dwarf masses estimated at 
the half-light radius with uncertainties om ■ The number of 
MW dwarfs in Walker et al. (2009) catalogue is iV obs = 23, 
where iV aa t is the number of satellite galaxies per realization 



with masses above mthres- This number strongly depends 
on the mass threshold, and for a MW-like halo it ranges 
between ~ 10 4 and 10 2 for m t hros = 10 7 and 1O 9 M , re- 
spectively. 

Interestingly, the masses of MW dSphs strongly con- 
strain the mass threshold for star formation in dark mat- 
ter subhaloes. We find that only those subhalos that were 
accreted with masses above 10 s and 10 9 Mq have average 
masses compatible with those in the sample of MW dwarfs 
with suitable kinematic data. Subhaloes that obey this con- 
dition may be identified with "luminous" satellites while the 
rest remain dark, in agreement with a broad number of theo- 
retical expectations (Bullock et al. 2009; Kravtsov 2010 and 
references therein). 

Comparing both panels we find that, independently of 
the host disc mass Md, satellite models with a cuspy inner 
profile provide a better match the MW dwarf galaxy masses 
than cored models, as expected from Fig. |11| 

However, there are a number of issues in this compar- 
ision that call for caution. For example, the census of MW 
dwarfs is still incomplete and limited by the surface bright- 
ness that observational surveys can reach. It is thus unclear 
the mass-size relation shown in Fig. [IT] holds for the entire 
satellite population. Also, we have implicitely assumed a lack 
of correlation between the mass and the stellar properties of 
dwarf galaxies. Although this is motivated by observational 
data (e.g. Mateo 1999; Walker et al. 2007; Penarrubia et al. 
2008a), self-consistent cosmological models that follow star 
formation/feedback in subhaloes are required to carry out 
a more quantitative analysis, as well as to explore in depth 
the contraints on galaxy formation models provided by the 
dSph population in the Local Group. 



5 SUMMARY 

We have used N-body simulations to explore the differ- 
ences that dark matter cusps and cores introduce in tidal 
mass stripping of satellite galaxies. Our models assume 
cosmologically-motivated initial conditions, where satellites 
are dark-matter dominated systems moving on eccentric or- 
bits within a much more massive host galaxy, which is com- 
posed of a "baryonic" disc embedded in an extended dark 
matter halo, whose density profile is assumed to be the same 
as that of satellite galaxies. Our findings can be summarized 
as follows 

• Our models show that the resilience of satellite galaxies 
to tides increases extraordinarily with the slope of the inner 
halo profile, 7. In particular, satellite models with 7 > 1.0 
cannot be fully disrupted by the tidal field of the parent 
galaxy, and always retain a bound remnant, even after losing 
more that 99.99% of mass to tides. In contrast, satellites 
with cored haloes 7 = 0.0 may undergo full tidal disruption 
if their orbits bring them in the vicinity of the disc. 

• Tides do not affect the inner structure of dark matter 
haloes, so that 7 stays constant during the evolution of satel- 
lite galaxies. Remarkably, the structural parameters of the 
satellite halo, e.g. the peak velocity u m ax and its location 
'"max, evolve in a manner that varies depending on the value 
of 7, but that is solely controlled by the total amount of 
mass lost to tides. 
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Figure 12. Mass threshold for star formation in satellite galax- 
ies, "i t hres! as a function of host disc mass. Points are grey-coded 
according to the accuracy of cuspy (upper panel) and cored (lower 
panel) satellite models in describing the mass-size relation shown 
by MW dwarfs (see Fig. |ll[ l. Note that (i) cuspy models provide a 
considerably better match to observations and (ii) a mass thresh- 
old at m t h res ~ 1O 8 -1O 9 M0 is clearly preferred for a large range 
of disc masses. 

To examine the differences that dark matter cores and 
cusps may introduce in the present satellite population sur- 
rounding spiral galaxies, we carry out ten different realiza- 
tions of merger trees that describe the hierarchical formation 
of a spiral galaxy with properties similar to those of the MW. 
These realizations are re-run adopting different values of 7 
as well as disc-to-halo mass ratios, Md/M v i r . We account 
for the fact that only a small fraction of subhaloes may be 
"luminous" by only considering those that fall into the host 
with masses above 10 s M@. Our main conclusions are the 
following 

• The net effect of tides is to decrease the average mass 
of satellite galaxies at all Galactocentric radii. Increasing 
the mass of the host disc clearly magnifies this effect: the 
average bound mass drops by a factor ~ 10 after injecting 
a disc component in the host with Md/M v i T = 0.1, and a 
factor ~ 30 for M d /M vil = 0.2 (see also D'Onghia et al. 
2009). 

• The present number of satellite galaxies that are ac- 
creted before re-ionization puts strong constraints on the 
inner profile of dark matter haloes. Our models predict that 
a MW-like galaxy may contain ~ 15 of these galaxies or- 
biting within the central 20 kpc if 7 = 1.0, in very good 
agreement with recent results from the Aquarius team (Gao 
et al. 2009). The number of these systems, however, quickly 
drops to zero in models that assume a cored halo profile 
and host disc component with Ma ~ 0.1M v ; r . These systems 
may be singled out by anomalous metallicity patterns, and 
thus represent an interesting observational case to unravel 
the inner structure of haloes. 



• In models that adopt a cored dark matter halo profile, 
discs efficiently deplete the inner regions of the host from 
satellites. In contrast, cuspy satellites can survive in the 
inner-most regions of spiral galaxies even after losing large 
fractions of their original mass to tides, which may have im- 
portant observational consequences for the detection of dark 
matter annihilation signals in gamma-ray surveys. 

• The relationship between half-light radius and enclosed 
dynamic mass shown by most Local Group dwarf galaxies, 
which spans approximately two orders of magnitude in size 
and three in mass, is successfully reproduced by our mod- 
els if adopting a 7 = 1.0 halo model. In contrast, cored 
halo models systematically underestimate the masses of the 
Ultra-Faint dSphs recently discovered in the MW. 

• In order to explain the size-mass relationship exhibited 
by the MW dSph population, a subhalo mass threshold for 
star formation is needed. We find that only those subhalos 
that were accreted with masses above 1O 8 -1O 9 M0 may have 
hosted luminous satellites that agree with the observational 
constraints, while the rest remain dark or beyond the present 
detection limits. 

• Our models show that high disc-to-halo mass ratios 
tend to lower the overall masses of satellite galaxies. A rela- 
tively massive M31 disc may explain why many of its dSphs 
with suitable kinematic data fall below the size-mass rela- 
tionship established from MW dwarfs. Recent estimates sug- 
gest that the disc of M31 may be a factor w 2 more massive 
than that of the MW (Hammer et al. 2007), lending support 
to this scenario. This explanation is likely to be validated (or 
challenged) by the results of many ongoing kinematic sur- 
veys of M31 dSphs, which will soon determine whether the 
difference in the size-mass relationship is systematic to the 
majority of the dSph population surrounding M31. 
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APPENDIX A: THE SEMI-ANALYTIC CODE 

To follow the dynamical evolution of dark matter subhaloes, 
we have constructed a simple semi-analytic code that com- 
putes their orbit, mass and structural evolution since they 
cross for the first time the virial radius of the host galaxy 
to the present. This code permits us to explore in a flexible 
manner aspects of galaxy formation that, if simulated with 
the aid of large N-body algorithms, would be notably expen- 
sive in terms of CPU requirements, such as the importance 
of cosmic variance in making observational predictions, or 
the effects of a varying the baryonic disc mass on the satel- 
lite galaxy population. 



Al Computation of orbits 

To calculate their orbits, we treat dark matter subhaloes 
as point masses that move in a host halo potential whose 
parameters evolve as shown in § |4.1| At each time step, our 
code solves the following equation of motion 



f = f h + i d + fd 



(Al) 



where the first term on the right side of the equation is the 
host halo force 



fh = -GM(< r)/r 2 - 
and 

rr 

M{< r) = 4tt / p(r')r' 2 dr' 



(A2) 



(A3) 



where p(r) has been defined in eq. [2] and [2]. The second 
term on the right side of the equation refers to the force 
induced by the baryonic disc, which we approximate as a 
Miyamoto-Nagai (1975) model (see 



4.11 



The third term on the right side of the equation is the 
dynamical friction force exerted by the diffuse sea of dark 
matter particles on massive satellites orbiting around their 
host. Several studies of satellite orbital decay have shown 
that, in spherical systems, Chandrasekhar's formula for dy- 
namical friction (Chandrasekhar 1943) is sufficiently accu- 
rate if the Coulomb logarithm is treated as a free parame- 
ter to fit to N-body orbits (e.g. Colpi, Mayer & Governato 
1999; van den Bosch et al. 1999). Semi-analytic methods 
that include Chandrasekhar's dynamical friction have been 
demonstrated to reproduce accurately the overall dynami- 
cal evolution of satellite galaxies in a cosmological context 
(e.g. Velazquez & White 1999; Taylor & Babul 2001; Zentner 
& Bullock 2003; Benson et al. 2004; Penarrubia & Benson 
2005) and therefore represent a useful tool for conducting 
extensive studies of a large parameter space. 

Dividing the background potential into a disc and dark 
matter halo components, Chandrasekhar's formula can be 
written as 



fdf — fdf.disc + fdf.halo — 
Vrol 



(A4) 



V Ie l) In A; 



' rcl 



where In Ad and In A^ are the Coulomb logarithms of the disc 
and halo components, respectively. The relative velocity of a 
satellite with respect to disc and dark matter halo particles 
is simply v rc i jd = v v d , , where = R\f d {Z = 0)| is 
the disc circular velocity on the galaxy plane [Z = 0); and 
Vrci.h = v, where v is the satellite velocity vector. 

For simplicity we assume here that the velocity distri- 
butions F(v) of both disc and host halo are Maxwellian and 
isotropic 



F(< w re i,i) = erf(Xi) 



(A5) 



where Xi = |vrei,i|/v2oi an d <J% is the one-dimensional 
velocity dispersion, which is defined as <7j(r) = 
l/ P i(r)Zpi(r')[f h (r') + f d (r')}dr'. 

Kazantzidis et al. (2004) have examined the inadequa- 
cies of the local Maxwellian approximation when applied 
to NFW halos and, although it may introduce problems to 
generate halos in perfect equilibrium, the deviations from 
Gaussianity of the exact distribution function have a negli- 
gible impact on the friction term (see their Fig. 4). 

Penarrubia et al. (2004) find that the best fit to self- 
consistent N-body orbits is obtained for a Coulomb loga- 
rithm In Ah = 2.1, in good agreement with the recent re- 
sults of Arena & Bertin (2007). Following Taylor & Babul 
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(2001), we use lnA^ = 0.5. We must emphasize, however, 
that these quantities are fairly sensitive to numerical as- 
pects, such as particle number and spatial resolution (e.g. 
Prugniel & Combes 1992, Wahde & Donner 1996) as well as 
to the host's density profile (Hashimoto et al. 2003, Just & 
Penarrubia 2005, Read et al. 2006, Arena & Bertin 2007). 
It is therefore necessary to note that the decay times of 
satellites orbiting in dark matter halos are uncertain to a 
significant degree, since orbital decay rates depend on the 
value assumed for Ad and A^. 

A2 Mass loss 

Satellite galaxies may lose a fraction of their original mass 
to tides if their orbits bring them into the inner regions of 
the parent galaxy. A simple way to estimate the bound mass 
fraction of a satellite galaxy is to define the tidal radius as 
the distance from the satellite centre at which the force due 
to the satellite's self-gravity and that from the external tidal 
force balance. For satellites moving on circular orbits, one 
can estimate that radius as (King 1962) 

^-( ^-d^ )^ (A6) 

where to is the angular frequency of the orbit and &h is the 
potential of the host halo. In calculating d 2 $h/dr 2 , we will 
average over the asphericity of the potential due to the disc 
component, so that 



d 2 $ h /dr 2 = d(- 



GM(< r)/r )/dr. 



(A7) 



If the satellite follows an eccentric orbit, one can still 
use eq. [A6| to calculate the instantaneous tidal radius at 
orbital pericentre, where tidal forces are strongest. 

A3 Evolution of the density profile 

As Fig. [5] illustrates, tidal mass stripping lowers the under- 
lying density of dark matter subhaloes at all radii, but more 
effectively in the outer-most regions. This process has been 
studied in some detail by Penarrubia et al. (2008b, 2009), 
who showed that the outer profile of dark matter haloes 
that lose a fraction (even as small as 10%) of its mass to 
tides can be well fitted by a power-law profile that scales as 
r~ 5 as r — > 00. In adopting the density profile outlined by 
eq. [2], this means that /3 — > 5. It is interesting to note that 
stripped haloes have f3 > 3, so that the resulting cumulative 
mass profile does not diverge at large radii. Remarkably, in 
the inner-most regions the dark matter halo profile is barely 
altered by the action of tides. In practical terms this means 
that 7 as well as a can be assumed to stay nearly constant 
throughout the dynamical evolution of our subhalo models. 

Furthermore, the evolution of the structural parameters 
of a subhalo as it is stripped can be accurately described by 
tidal tracks, i.e. mono-parametric functions that only depend 
on the fraction of mass lost to tides (see Fig. |6|. 

We can use the above information to describe the 
changes induced by tidal mass stripping on the density pro- 
files of individual dark matter subhaloes as a function of 
their bound mass fraction. This is done in a fairly simple 
way: (i) for m s /m s (t — 0) < 0.9 the indices of eq. [2] are 
modified as (a, ,8,7) — > (0,5,7). (ii) F° r a given value of 
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Figure Al. Evolution of the density profile of our satellite 
galaxy model moving on a highly eccentric orbit (r per i : r apo = 
1 : 100) at different orbital apocentres (Fig. |5j. Profiles are ar- 
bitrarily normalized by the quantity p n orm, which is kept fixed 
in each panel to emphasize the density drop resulting from tidal 
mass stripping. Dotted lines correspond to our analytical model 
predictions (see text). 



m a /m s (t = 0) we calculate (r^ax) u^ax) using eq. [Ij. (iii) 
The new scale radius r' s is then calculated from the condi- 
tion dV 2 /dr| r=r j nax = 0, (iv) whereas p'o is obtained using 
the definition 14(r^ ax ) = t4 ax . 

This simple method permits us to follow the evolution 
of the density profile of satellite galaxies as they orbit within 
their host with an appropiate degree of confidence. For ex- 
ample, symbols in Fig. |A1| illustrate the density profiles 
shown in Fig. [5] whereas dotted lines show the profiles de- 
rived from the above recipes. Except for the excess of mass 
visible at large radii in some cored models, which corre- 
sponds to unbound material on its way to escape from the 
subhalo, our semi-analytic algorithm successfully matches 
the evolved density profiles, even those of subhaloes that 
have lost a large fraction of their initial mass. This code will 
permit us to explore in a very efficient way aspects of galaxy 
formation that would otherwise require an extensive use of 
supercomputers, such as accounting for cosmic variance to 
make theoretical predictions, or the effects that varying the 
host galaxy disc-to-halo mass ratio would induce on the sur- 
viving satellite galaxy population. 



A4 Tests 

We have used the N-body simulations of ^3] to test our mass 
loss algorithm, as well as the analytical description of the 
effects that tides cause on the density profiles of dark mat- 
ter haloes. The results are shown in Fig. |A2| where we plot 
the mass our subhalo models as a function of time as de- 
rived from our N-body models (see Fig. [3]l , as well as from 
the above semi-analytical method (dotted lines). The agree- 
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Figure A2. Mass evolution of our satellite galaxy models mov- 
ing on different orbits. Dotted lines correspond to our analytical 
model predictions (see text). 



merit between both techniques is reasonable, taking into ac- 
count the extreme mass-loss events considered here and the 
large range of orbits and structural parameters adopted in 
this study. More importantly, our simple semi-analytical de- 
scription is able to capture the strong effect that discs have 
on the mass evolution of dark matter subhaloes depending 
on the inner shape of the density profile, which will per- 
mit us to explore the effects that baryonic discs have on the 
satellite galaxy population depending on whether haloes are 
cuspy or cored. 
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